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ABSTRACT 

We argue that inner rings in barred spiral galaxies are associated with specific 2D and 
3D families of periodic orbits located just beyond the end of the bar. These are families 
located between the inner radial ultraharmonic 4:1 resonance and corotation. They 
are found in the upper part of a type-2 gap of the xl characteristic, and can account 
for the observed ring morphologies without any help from families of the xl-tree. Due 
to the evolution of the stability of all these families, the ring shapes that are favored 
are mainly ovals, as well as polygons with 'corners' on the minor axis, on the sides of 
the bar. On the other hand pentagonal rings, or rings of the NGC 7020 type hexagon, 
should be less probable. The orbits that make the rings belong in their vast majority 
to 3D families of periodic orbits and orbits trapped around them. 
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1 INTRODUCTION 

Rings are often observed in barred galaxies and can be 
spectacular structures. They come in three kinds: nuclear 
rings, near the nucleus of the galaxy, inner rings, surround- 
ing the bar, and outer rings, of a considerably larger di- 
ameter. In this pap er we focus our attention to the study 
of the inner rings. iButal lll995l . hereafter B95) has made 
a statistical study of 3692 ring galaxies in the southern 
hemisphere, as well as more in depth studies of specific 
objects (e.g. iButa et alJl200ll Il999t iButa fc Pureed Il998t 



iButa. Purcell fc Cro cker 199fi|). He found that inner rings 
are a frequently encountered feature of barred galaxies. 
Their mean axial ratio is 0.84 ± 0.10, elongated along the 
bar major axis. This ratio, however, varies from one galactic 
type to another as w ell as among different ring morpholo- 
gies (see Table VII in lButa fc Combeslll996ri . In some cases 
of typical inner rings the axis ratio can be as low a s 0.63 
(NGC 6782) or even 0.49 (A0106.7-3733) teutalll986l) . B95 
mentions a dozen galaxies for which there could be an in- 
trinsic misalignment between the bar and the inner ring. Be- 
cause of projection effects, however, spectroscopic measure- 
ments are necessary in order to confirm this misalignment. 
In the vast majority of cases the inner ring encircles the bar 
and touches its extremities. There are, however, cases where 
the ba r underfills what appears to be an inner ring. IButal 
( 1986) mentions 30 such cases in his sample of 1200 objects. 
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Two typical examples are also given by B95, namely NGC 
7098 and NGC 3450. 

The link between rings and resonances has be en made 
already in the '70s (e.g. lSchommer fc Sullivanll976ll . It took, 
however, the advent of test particle hydrodynami c simula- 
tions i n the ' 80s to establish it fully. In particular. ISchwarzl 
dl98llll984ah followed the response of gas, modeled by sticky 
particles, to a bar forcing. He found that ring structures 
can form, and that, as in real galaxies, they can be classi- 
fied into outer, inner and nuclear. Schwarz linked the outer 
rings to the outer Lindblad resonance (OLR), the inner 
ones to the inner 4:1 resonance or inner ultraharmonic res- 
onance (iUHR), and the nuclear ones to the inner Lind- 
blad resonance (ILR). The first statistical arguments came 
from lAt hanassoula et a ll <ll982l) who, using the sample of 
Ide Vaucouleurs fc Butal Jl98(J) . showed that the ratios of in- 
ner to outer ring sizes is compatible with the outer ring be- 
ing at the OLR and the inner one at corotation, or at the 
iUHR. This statistics has been since repeated (B95), using 
much better observational samples, but the result did not 
change. 

Rings form by gas accumulation at resonances (e.g. 
lButalll999l . and references therein), and this explains their 
predominance on images of galaxies which show best the 
sites of population I objects. It is, however, now clear 
jAthanassoulalll992bl) that the gas response is determined 
to a large degree by the periodic orbits in the underlying 
gravitational potential. For this reason, understanding the 
orbital dynamics of rings is essential not only for under- 
standing their morphology per se, but also for understanding 
properties related to the gas component, e.g. star formation. 
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One of the main uses of orbital structure studies is that 
they provide information on the orbits that are the backbone 
of the various galactic structures, and thus on their mor- 
phologies. Studying the orbital structure in an appropriate 
Hamiltonian system we get the periodic orbits that could be 
responsible for the appearance of morphological features of 
real galaxies. We are thus interested only in stable periodic 
orbits, since they trap around them the regular orbits of the 
system. While in 2D bar models the ba ckbone of the bar i s 
a single family, namely the family xl iContopoulosi ri983). 
in 3D models the orbits supporting the bar are related to 
a whole t ree of families of pe riodic orbits, which we called 
'xl-tree' iSkokos et al . 2002a). Under the assumption that 
bars are fast, their end should be close and within corota- 
tion <|ContoDouloslll980h . Due to the fact that inner rings 
surround the bar, one intuitively links them to stable orbits 
not belonging to the xl-tree, but occupying the area just 
beyond the bar in the real space. 

Although an association of a resonance location to an 
inner ring can be easily done in a first approximation, there 
is so far no detailed study relating the various observed mor- 
phologies of inner rings to the orbital dynamics at their re- 
gion. As Buta remarks "...ring s include in formation on ... the 
properties of periodic orbits" feutall2002l) . This information 
has not been registered until now, especially for 3D barred 
models, and this is exactly what gave us the motivation for 
the present paper. 

Inner rings can be either smooth, in which case we refer 
to them with the generic name of 'oval', or they can have 
characteristic 'angles', or 'kinks', or 'corners' in their outline. 
These can be more or less strong and, depending on their 
strength, may give the ring a polygonal, rather than oval, 
geometry. They occur at characteristic locations along the 
ring. Thus such 'corners', or 'density enhancements' often 
occur near the bar major axis, giving the ring a 'lemon-like' 
outline. In other cases, we observe that the ring is squeezed 
along the major axis of the bar, so that the tips of the 
'lemon' structure vanish, the ring becomes rather polygonal- 
like with sides roughly parallel to the bar's minor axis at its 
apocentra, and 'corners' cl ose to the minor axis. T ypical ex- 
amples are UGC 12646 fe.g lButa fc Combesll996l - Fig. 17) 
IC 42 90 (see the distribution of the HII regions in | Buta et all 
I199SI. - Fig. 8), NGC 3351 (e.g. ISandaee fc Bedkdll994 
panels 168, 170). In such cases it is reasonable to speak 
roughly about a hexagon with two sides parallel to the minor 
axis of the bar. The above mentioned objects (UGC 12646, 
IC 4290 and NGC 3351) are typical examples. Nevertheless, 
in most cases all these features characterize only parts of 
the rings and perfect symmetry, e.g. with respect to the bar 
major axis, is not always observed. 

Despite the complexity of the observed structures of the 
inner rings, and the fact that they are located outside the 
bar, they have been in general vaguely associated with or- 
bits at the inner 4:1 resonance <ISchwarzlll984bl: lButalll99l 
2002 ; ). ISchwarzl Jl984b) had already realized the presence 
of squeezed ovals with corners. He was, however, able to 
reproduce them only in the c ase of weak b a rs, or by in- 
voking a lens-like component <ISchwarzlll985r ). iButal (2002, 
- Fig. 3) invokes combinations of diamond- with barrel-like 
orbits to explain ring morphology. In particular barrel-like 
orbits have been considered necessary in order to provide 
cloud-cloud collisions on the sides of the bars (the regions 



where the cloud-cloud collisions happen are indicated with 
"C" in ISchwarzlll984bl) . These orbits help in the formation 
of rings which have a hexagonal rather than diamond-like 
geometry and have two of the sides more or less parallel to 
the minor axis of the bar. 

Note that neither every diamond-like, nor all hexago- 
nal orbits are necessarily related with inn er rings. E.g. the 
planar diamond-lik e xl orbits in Fig. 62 inlButa fc Combed 
i|l996li , taken from IContopoulos fc Grosb0ll " i^89l) ~ are far 
away from the end of the bar, which is close to corotation. 
We also note that hexagonal rings with cusps on the major 
axis of the bar, and thus with two sides parallel to the major 
axis of the bar, are rare. There is only a notab le example o f 
this morphology, namely NGC 7020, studied bv lButal lll99Cl) . 

In the present paper we study the orbital behavior at 
the region where the appearance of inner rings is favored. 
We do this in the case of 3D Ferrers bars. Studying the en- 
ergy width over which the various families of periodic orbits 
exist, their morphology, and mainly their stability, we iden- 
tify the orbital behavior in four archetypical morphologies 
of inner rings. In section|5]we present these structures in the 
cases of four galaxies. In section0we give a brief description 
of the model we use for our orbital calculations, while the 
families which build the rings are described in section ^] Fi- 
nally we discuss our results in section 03 and enumerate our 
conclusions in section [S] 



2 RINGS MORPHOLOGY 

Let us first introduce four inner ring morphologies. Two of 
them are typical ovals and the two other are exceptional 
cases of inner rings. Although the images 1 shown in Fig. 
are not of high resolution, they are nevertheless able to 
demonstrate the features we refer to. Most inner rings are 
oval, sometimes with a somewhat lemon shape because of 
density enhancements at the bar major axis. We illustrate 
this latter type in Fig. with an image of NGC 6782. The 
isophotes peak on the major axis, indicated with arrows on 
the second outermost isophote, giving to the ring a cuspy 
or pointed oval morphology. There are no "corners" close to 
the minor axis of the bar, which would have given to the 
ring a diamond shape. In Fig. Et> we show IC 4290 which 
has an oval ring with characteristic breaks or corners. These 
are indicated with arrows on the overplotted isophotes. The 
outermost isophote is again oval-shaped and even slightly 
pointed close to the major axis. The oval, however, is not 
smooth and forms, or gives the impression of, angles. One 
can also discern a structure in the form of an "I" with the 
vertical part of the symbol being the bar. The two other rings 
given in Fig.^are exceptional cases. The ring in Fig. [Hi has 
an almost pentagonal geometry. It is observed in an image 
of NGC 3367. It can be described as a pseudo-ring, since the 
ring appears to be made of segments of spiral arms. F inally 
in Fig. 01 we give a deprojected image of NGC 7020 ( Buta 

1 The image of NGC 6782 originates from an ScTScI Hubble 
Space Telescope image. It can be viewed in high resolution at 
several web sites. The three others are DSS images. The image 
of NGC 7020 in Fig. has been deprojected. The deprojected, 
as we l l as o ther images of the NGC 7020 hexagon can be seen in 
IButal Il99(f> 
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Figure 1. Examples of the four ring morphologies discussed in 
the paper, (a) The oval with cusps on the major axis (NGC 6782), 
(b) a typical oval-polygonal ring with sides parallel to the minor 
axis of the bar close to the ends of the bar (IC 4290), (c) a pen- 
tagonal ring in NGC 3367 and (d) the hexagonal ring with cusps 
on the major axis (NGC 7020). At the bottom left corner of each 
panel there is a sketch indicating schematically the corresponding 
morphology. 



1990), which has a type of hexagonal ring with cusps on the 
major axis of the bar and two sides parallel to it. The "cor- 
ners" are indicated with arrows. Of course the geometrical 
shapes we use above are idealized and represent the com- 
plicated morphology of inner rings only schematically. They 
are, nevertheless, most useful as guiding lines, in the same 
way as perfect ellipsoids are useful for describing bars. 



3 MODEL 

To calculate the 3D ring or bits we use the fiducia l case of the 
model described in detail in ISkokos et aP(l2002al) . It consists 
of a Miyamoto disc, a Plummer bulge and a Ferrers bar. T he 
potential of the Miyamoto disc jMivamoto fc Nagaill975T) is 
given by the formula: 



GM D 



y/x 2 + y 2 + {A + /B2~kP 



(1) 



Mb represents the total mass of the disc, G is the gravita- 
tional constant, and A and B are scalelengths such that the 
ratio B/A gives a measure of the flatness of the model. 

The bulge is a Plummer sphere, i.e. its potential is given 

by: 



$s = 



GM S 



yjx 2 + y 2 + z 2 + e 2 



(2) 



where e s is the bulge scale length and Ms is its total mass. 

Finally, the bar is a triaxial Ferrers bar with density 
p{x): 



p{m) = 



W5Mb 2,2 , . , 

^-^(1-m) for m<l 



for m > 1 



where 

2 y 2 x 2 z 2 

m = ^t + tt + ^7, a> b> c. 

a z tr cr 



(3) 



(4) 



In the above a, b, c are the principal semi-axes, and Mb is 
the mass of the bar component. The correspondi ng potential 
$b and the forces are given in IPfennigerl lll984h in a closed 
form, well suited for numerical treatment. For the Miyamoto 
disc we use A=3 and B=l, and for the axes of the Ferrers bar 
we set a:b:c = 6:1.5:0.6. The masses of the three components 
satisfy G(M D + Ms + M B ) = 1. We have GM D = 0.82, 
GM S = 0.08, GM B = 0.10 and e s = 0.4. 

The length unit is taken as 1 kpc, the time unit as 
1 Myr and the mass unit as 2 x 10 M®. The bar rotates 
with a pattern speed f2{,=0.054 around the z-axis, which 
corresponds to 54 km sec" 1 kpc -1 , and places corotation at 
6.13 kpc. 



4 ORBITS 

The Hamiltonian governing the motion of a test-particle in 
our system can be written in the form: 

H = ^(pl+pl+pl) + V(x,y,z) - Q b (xp y - yp x ), (5) 

where p x , p y , and p z are the canonically conjugate momenta 
of x, y and z respectively and V(x, y, z) is the total potential 
of the combined three components of the model: disc, bar 
and bulge. We will hereafter denote the numerical value of 
the Hamiltonian by Ej and refer to it as the Jacobi constant 
or, more loosely, as the 'energy'. 

A periodic orbit in a 3D Hamiltonian system can be ei- 
ther stable (S) or exhibit various types of instability. It can 
be simple unstable (U), doubl e unstable (D) or complex 
unstable (A) (for definitions see IContopoulos) 2002, chapter 
2.11). The method used for findin g periodic orb i ts and their 
stability is described in detail in ISkokos et aP ll2002al) and 
in references therein. Periodic orbits are grouped in fami- 
lies, along which the energy Ej changes and the stability of 
the orbits may also change. The S — > U transitions of peri- 
odic orbits are of special importance for the dynamics of a 
system, since in this case a new stable family is generated 
by bifurcation. A useful diagram for presenting the various 
periodic orbits of t he system is the 'characteristic' diagram 
<Contopoulojl200l section 2.4.3). It gives the x coordinate 
of the initial conditions of the periodic orbits of a family as 



a function of their Jacobi constant Ej . In the case of orbits 
lying on the equatorial plane (z — p z = 0) and starting per- 
pendicular to the x axis (p x = 0), with y = and p y > 0, we 
need only one initial condition, x, in order to specify a pe- 
riodic orbit on the characteristic diagram. Thus, the initial 
conditions of such orbits are fully determined by a point on 
the characteristic diagram. Initial conditions of orbits which 
stay on the equatorial plane but do not start perpendicular 
to the x axis and of orbits which do not stay on the equa- 
torial plane are not fully defined by their x coordinate. In 
these, more general cases, the representation on the (Ej,x) 
diagram is not very useful. 
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In our model, the highest energy orbits that support 
the bar can be found at the region of the rad ial 4:1 res- 
onance ijSkokos et al]l2002al IPatsis et al.ll 20031) . Since the 
rings surround the bars, we first look for families of periodic 
orbits beyond the energy which corresponds to this reso- 
nanc e and we follow their stability as we approach corota- 
tion. lAthanassoulal lll992al) has found several families of or- 
bits beyond the 4:1 resonance in 2D Ferrers bars (see Figs. 2 
and 3 of that paper). We did not know, however, which of 
those could trap matter around them and support ring struc- 
tures, because their stability was not studied in that paper. 
We address this problem here in the more general case of a 
3D Ferrers bar. 




4-0.1 The f- group 

In all 3D Ferrers bars studied bv lSkokos et al 1 J2002allbh the 

characteristic of the xl family has x values which increase 
with Ej until a maximum is reached, at which point the x 
values start decreasing. Beyond the local maximum we find 
the characteristics of other families, which for larger Ej have 
initial x's which are l arger than those of xl. In othe r words, 
we have type-2 gaps JContopoulos fc Grosbollll989lb In 2D 
models of barred galaxies the f amily found beyond the local 
maximum is the family xi(2) l)ContoDoulosll2002T Fig. 3.8b 
in). In our 3D models the first family we find beyond the 
gap is a 2D 4:1 family, which we call f. Its morphology is 
rhomboidal as the mor phologies of the corresponding fami- 
lies of the 2D models in lContopoulos fc Grosbol <ll989l) and 
lAthanassoulal (ll992aT ) . However, in the 3D model, it is ac- 
companied by a forest of bifurcating families. We call this 
group of families of periodic orbits the 'f-group'. 

In Fig. |21 we give the (Ej,x) characteristic of f, which is 
the small curve above the xl characteristic indicated with f. 
The lower part of the curve, i.e. smaller x for the same Ej, is 
composed mainly of stable f orbits. Nevertheless, instability 
strips do exist, but are not important. The upper part of the 
curve is mainly composed of unstable orbits. In the figure 
black segments indicate stability and grey instability. For the 
sake of clarity we do not include in the figure the projections 
of the characteristics of the 2D and 3D families bifurcating 
from f. 

The 'mother' planar 2D f-family starts existing at about 
Ej ~ —0.20166, and the morphology of its orbits is as in 
Fig-EK- It remains basically stable until Ej ~ —0.198, with 
only small instability regions in between. With increasing 
energy, the morphology of the f orbits changes and becomes 
hexagonal (Fig.|3)). We thus have a smooth transition from 
a basically 4:1 to a 6:1 morphology. The two points marked 
"A" and "B" on the stable part of the characteristic indi- 
cate the location of the two f-orbits shown in Fig. |2K and 
Fig-Et) respectively. Beyond Ej ~ —0.198 the orbital multi- 
plicity, i.e. the number of intersections of the periodic orbit 
with the x-axis before it closes, changes. At an S— >U-^S 
transition close to Ej ~ —0.2004 we have two 5:l-type 2D 
families bifurcated from f. Namely, frl (Fig. Efj) at an S— >U 
transition, and fr2 (Fig. EJi) at the nearby U— >S transition. 
We note that frl orbits do not start perpendicular to the 
x axis. For these and all the rest of the bifurcating families 
in the paper w e follow the nomenclature rules introduced 
in lSkokos et al . ( 2002a). Soon after its bifurcation, however, 
the stable frl family becomes unstable, while the initially 



Figure 2. Characteristic diagram of families xl, f and s. We 
observe a typical type-2 gap at the radial 4:1 resonance region. 
Stable parts are indicated with black lines and unstable ones with 
grey lines. The capital letters A to F, as well as U, indicate the 
location of specific orbits which are discussed further in the text. 
At this diagram L4 5 would be at rj (—0.1955,6.0). 




Figure 3. 2D f-orbits. In (a) and (b) we have stable f orbits 
located at points "A" and "B" on the f characteristic in Fig. [5] In 
panels (c) and (d) we give typical stable frl and fr2 orbits. In this 
and subsequent figures the bar major axis lies along the j/-axis 
and the units are in kpc. 



unstable fr2 becomes stable. The f orbits along the unstable 
part of the characteristic are diamonds with cusps on the 
axes (Fig. [IJ. These orbits are characterized by segments, 
which resemble straight line sides. 

Stable 3D families are introduced at S^U transitions 
of family f at Ej « -0.2016, -0.2002 and -0.1993, and are 
the families fvl, fv3 and fv5, respectively. These families, to- 
gether with their stable bifurcations, contribute significant 
stable parts to the area. The three projections of five or- 
bits of 3D families of the f-group are given in Fig. |S] The 
two orbits at the top of the figure belong to family fvl, for 
Ej = —0.201 and —0.181 respectively, and show that the 
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-6 -4 -2 2 4 6 

Figure 4. A 2D unstable f-orbit from the upper part of the f- 
characteristic (drawn with a grey line in Fig. |5J . 



morphological evolution of the (x, y) projection of this fam- 
ily follows the morphological evolution of the 2D family f. 
The third and fourth rows from the top give the morphol- 
ogy of orbits from two bifurcations of family fvl, which in 
turn has bifurcated from family f. They are introduced in 
the system as stable, but do not have large stability parts, 
since they become soon unstable. The fifth and the sixth 
orbits from the top belong to family fv3, and are given for 
Ej = —0.200 and —0.193, respectively. The morphology of 
their (x, y) projections is again similar to the morphology of 
a 2D f-type family (Fig.|^K,b) of the same energy. Finally the 
last orbit is an example of an orbit from family fv5, close to 
its bifurcating point. We note that fv5 becomes complex un- 
stable at Ej — —0.19697. At this energy it has already large 
deviations above and below the equatorial plane, i.e. large 
| z | values, and thus it would not contribute significantly to 
the local disc density, even if it was stable. 



4-0.2 The s-group 

In Fig. |2]we see that, below and to the right of the charac- 
teristic of f, there is another, longer curve, roughly speaking 
parallel to the decreasing branch of xl. This is the charac- 
teristic of a 2D 6:1 family, which we call "s". It exists for 
Ej > —0.1990. As in the case of the f-group, we have here 
a group of families of periodic orbits (s-group). The group 
includes s and several 2D and 3D families bifurcated from 
it. The stable part is the one with the larger x's for the same 
energy. As in the previous case, along the stable part we have 
a morphological evolution of the orbits and tiny instability 
strips, where new families are born. We can follow the mor- 
phological evolution of the s family of orbits in Fig. |S] from 
(a) to (b) and to (c). They correspond to the points labeled 
on the characteristic of s in Fig. [5] with the capital letters 
"C" , "D" and "E" , respectively. Along the characteristic of 
s we have a smooth transition from a basically 6:1 to an 8:1 
morphology. An S^U transition at Ej « —0.19876 intro- 
duces in the system the planar stable family srl (Fig. HUl), 
which is a 7:1 type orbit. Along the long unstable part of 
the s characteristic, orbits are hexagonal with cusps on the 




! v I 



-2 2 4 6 -6 



■2 2 4 E ■£ 



-3 2 4 e 




XJJ 



1 V 1 



■6 -4- -2 2 4. fl -6 -A -£ £ 4, 4 



4 4 3 4 6 




fvl.l 



-B-4-E02 4S£-4 £ QJ* 5-6-4-2 024 




.l.l.l 



6 -4 -2 P 3 1 B -6 -1 -2 £ 4 



i T. 2 4 



W 3 



■S5S4 6-«-4-2fl24fl 




(XKXD0- 



ft 3 



-4 -f & 2 4 



-4-2 i 




Fv5 



.2 a 24&-S-<.2iD 2 4 6 6 



Figure 5. Stable 3D f-orbits. The names of the families are noted 
on the right side of each set of projections. The projection planes 
are indicated in the upper left corner of each frame. 



major axis of the bar (Fig. HjJ:). Family s has another sta- 
ble part at the largest Ej's and small x values (see Fig. 
However, the morphology of the s orbits at these energies is 
characterized by big loops (Fig. 0), so that they cannot be 
used to build rings. 

At three other S^U transitions with tiny unstable 



parts, at Ej 



-0.19892, ir- 



respectively, three stable 3D families bifurcate. From smaller 
to larger energies, we call them svl, sv3 and sv5. Fig.|7|gives 
the morphology of these families. In Fig. |7Ji we observe the 
three projections of svl close to its bifurcating point. In 
Fig. [TJj we give sv3 also close to its bifurcating point and 
in Fig. [TJj at larger energies. Finally Fig. |7|I shows sv5 just 
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Figure 6. 2D orbits of the s-group. From (a) to (b) and then to (c) we see the evolution of the orbital morphology along the stable 
branch of the s characteristic as we move from "C" to "D" to "E" (Fig.|5J- In (d) we observe a planar stable srl orbit, in (e) an unstable 
s orbit (point U in Fig. [2), and in (f) a stable s orbit with loops (point F in Fig. |2*1, 



after its bifurcation from s. The morphological evolution of 
the face-on projections of the 3D families resembles the evo- 
lution of the stable s orbits as the energy increases. Both sv3 
and sv5 families have complex unstable parts, starting from 
E 3 « -0.19424 and Ej ^ -0.1972 respectively. However, 
as in the case of fv5, this happens at energies where \z\ is 
large, and thus the orbits of these families at these energies 
are not significant for the local disc density. 



5 DISCUSSION 

In 3D Ferrers bars, families of periodic orbits beyond the gap 
at the 4:1 radial resonance, where the characteristic of the 
planar xl family has a local maximum (Fig. I^J^Jire intro- 
duced through tangent bifurcations (see IContopo ulos 2002, 
pg.102 for a definition). Tangent bifurcations are fundamen- 
tal to the study of nonlinear systems since they are one of 
the most basic processes by which families of periodic orbits 
are created. In 2D barred potential s they correspond to the 
gaps at the even radial resonances llContopoulos fc Grosboll 
1989). In such a bifurcation a pair of families of periodic 
orbits is created "out of nothing" . They are also referred to 
as saddle-node bifurcations. One of the newborn sequences 
of orbits is unstable (the saddle), while the other is stable 
(the node). A characteristic diagram of this type of bifur- 
cations is shown schematically in Fig.[S] The characteristics 
of families f and s in Fig. are of this type, which means 
that these families are not connected to the xl-tree. Further- 
more, they build their own group of families, i.e. their own 



trees. Despite the fact that the stability is not indicated, one 
can see from the sha pe of the characteristic s that in the 2D 
Ferrers bar model in lAthanassoulal <ll992aT) the families in 
the upper part of the 4:1 resonance gap are introduced via 
tangent bifurcations. 

Close to the Ej minimum at which the f and s families 
appear, the orbits are characterized by 4:1 and 6:1 mor- 
phologies, respectively. As we move along the stable branch 
towards higher values of Ej, the f orbits initially become 
rounder, and then for yet larger Ej , they develop cusps along 
the minor axis of the bar, evolving morphologically, roughly 
speaking, to a hexagon. Note that this hexagon has two sides 
parallel to the minor axis of the bar. The morphological evo- 
lution of the stable s orbits forms a similar sequence, now 
from a 6:1 to an 8:1 morphology, again via an ovalish shape. 
Note, however, that the geometrical shape of the orbit de- 
picted in Fig.|SJi changes very fast, i.e. the s orbits have this 
morphology only in a very tiny energy interval. It evolves to 
a lens-type morphology and finally to the octagonal shape 
we see in Fig. [(J, which develops loops along the major axis 
of the bar. In order to visualize a collective effect of the pres- 
ence of such octagonal orbits in the system, we plot in Fig.|U] 
successive orbits of this type selected at equally spaced in- 
tervals in energy. This could lead to an octagonal ring or, 
most likely, to a hexagonal ring with additional density en- 
hancements close to the major axis of the bar. The shapes of 
the rings that are favored due to the stability of the planar f 
and s orbits can be described in general as ovals. Hexagonal 
morphologies support motion parallel to the minor axis at 
the apocentra of the rings. 



On the 3D dynamics and morphology of inner rings 7 




5450246 64202 + 6 6420246 




6420246 64202 + 6 6420246 




S 4 2 2 4 6 64202 + 6 6420246 




-6 -4 -2 2 4 G -6 -4-2 2 + 6 -6-4 -2 2 4 6 



Figure 7. Stable 3D orbits of the s-group. (a) orbit from the svl 
family, (b) sv3 orbit close to the bifurcating point, (c) sv3 orbit 
away of the bifurcating point, (d) orbit from the sv5 family. 




Figure 8. Schematic representation of a tangent bifurcation. The 
family starts existing for energies larger than E m i n . 
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Figure 9. A set of stable, octagonal s orbits. Note that, while 
the outer boundary is octagonal, the inner boundary tends to a 
hexagonal shape. 



The 2D radial bifurcations from f and s are associated 
with the 5:1 and 7:1 resonances, respectively. The corre- 
sponding 5:1 and 7:1 families come in pairs, with orbits sym- 
metric with respect to the bar major axis. Thus, although 
they do have stable parts, they do not in general support a 
particular morphology if we consider both of their branches. 
Indeed, two 5:1 symmetric periodic orbits together represent 
a 10:1 morphology. In this case, non-periodic orbits trapped 
around them will practically contribute to the formation of 
a smooth ring around the bar. In cases where only one of 
the two branches is populated, the ring will have a mainly 
pentagonal or heptagonal geometry. 

The existence of the 3D families is very important for 
the prevalence of a ring structure around the bar, because 
they extend the volume of the phase space occupied by sta- 
ble ring-supporting orbits. Such families are introduced at 
the S— >U transitions of our model and some of their mem- 
bers are plotted in Figs. |S] and |7| Their vertical thickness in 
the z direction remains low for considerable energy intervals, 
at least for the orbits with face-on projections without cusps 
or loops (Fig. . This is in good agreement with observa- 
tions, since inner rings have never been seen to stick out 
of the equatorial plane in galaxies which are observed edge- 
on. These new families remain practically stable for energy 
intervals over which their face-on projections, and the cor- 
responding f or s 'mother' families support the ring. On the 
other hand, over the energy intervals where these families 
are complex unstable, their face-on projections have loops 
and their average height above the equatorial plane has also 
increased. Since they are complex unstable they are not as- 
sociated with further bifurcating families (for definitions re- 
lated with instabilities in 3D systems see e.g. IContopoulosI 
2002). The presence of complex instability introduces chaos 
in the system, however at energies for which the families 
have members reaching large distances above and below the 
equatorial plane. Thus their importance for the density of 
the rings is already small. 
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The two groups of f and s families of periodic orbits 
include practically all relevant morphologies of orbits one 
can find between the radial 4:1 resonance and corotation. 
Although we can not exclude the existence of n:l (n>8) 
type families, we didn't find any such orbit remaining close 
to the equatorial plane. Since, furthermore, the f and s family 
orbits are sufficient on their own to account for the rings, 
we do not discuss the n:l (n>8) orbits any further. 

5.1 Interpretation of the ring morphologies 

In this section we will show that inner ring morphologies 
can be built using combinations of orbits of the stable 2D 
and 3D families which we described in the previous sections. 
Let us start with the most common morphologies. There is a 
large number of orbital combinations allowing the formation 
of the types of rings illustrated in Figs. ^ and b. The mor- 
phology of the s orbits that have an outline as that shown 
in Fig. |U renders well that of the most frequently encoun- 
tered inner rings (Fig. EJi). This is still true if we add to 
them stable svl and sv3 orbits, as those depicted in Figs.|7^ 
and b, respectively. In order to demonstrate this, we give in 
Fig. 110b the face-on profile of the weighted oval- or lemon- 
shaped orbits of families s for —0.1990 < Ej < —0.1985, svl 
for -0.1989 < Ej < -0.1968 and sv3 for -0.1987 < Ej < 
-0.1982 and -0.1976 < Ej < -0.1958 (sv3 is simple unsta- 
ble for -0.1982 < Ej < -0.1976). These intervals have been 
chosen so as to exclude all unstable orbits and all orbits with 
\z\ > 0.8. We consider orbits with \z\ > 0.8 as contributing 
little to the rings density. We use the s ame t echnique as in 
iPatsis. Skokos fc Athanassoulal \ 2001 1 2003ft . Namely, we 
first calculate the set of periodic orbits which we intend to 
use in order to build the profile. Then we pick points along 
each orbit at equal time step s. The 'mean density' of each 
orbit (see lPatsis et alJl 2002L - section 2.2) is considered to 
be a first approximation of the importance of the orbit and is 
used to weight it. We construct an image (normalized over 
its total intensity) for each calculated and weighted orbit, 
and then, by combining sets of such orbits, we construct 
the weighted profile. The selected stable orbits are equally 
spaced in their mean radius. The step in mean radius is the 
same for all families in a figure. 

In Fig. 110b we observe that all orbits of the s, svl and 
sv3 families are confined inside a very narrow ring. Fami- 
lies svl and sv3 are 3D. The orbits considered here are also 
vertically confined to a thin layer of a few hundreds pc hav- 
ing in the middle the equatorial plane of the galaxy model. 
For this, and all subsequent similar figures, we give in TabQ 
the energy intervals from which we have taken the orbits of 
a family in order to construct the profile. For building the 
profile we used orbits from all the available energy intervals. 
Due to the orbital crowding inside the ring, the orbits of 
the three families intersect each other, and thus the corre- 
sponding gas flow will be characterized by the presence of 
numerous cloud collisions in the same area. 

If we add, on top of this, orbits trapped around f, fvl 
and fv3, at the appropriate energy intervals, the morphol- 
ogy of the ring does not change significantly. In Fig. 110b we 
add stable orbits belonging to the families f for small ener- 
gies (-0.2016 < Ej < -0.1993), fvl for -0.2016 < Ej < 
-0.1974 and fv3 for -0.2002 < Ej < -0.1967. Again orbits 
with \z\ > 0.8 are excluded. The total energy width over 



which the oval-building orbits exist and specify the mor- 
phology of the inner ring is AEj — 0.0058. The contribution 
of 2D families is over a width AEj = 0.0028, while the 3D 
families support rings in an interval AEj — 0.0058. This in- 
dicates that 3D orbits are essential in building inner rings. 
It is also evident that the width of the ring in the face-on 
view is set by the contribution of orbits belonging to the f- 
tree. These orbits could be the origin of broader inner rings 
observed in the near-infrared. 

In Fig. HOt we apply a gaussian filter to the image of 
Fig. 110b . in order to show, in a first approximation, the 
shapes of secondary features which could be supported by 
the orbits discussed here. The colour bar below it shows the 
correspondence of colours and surface density on the ring. 
As we move to the right on the colour bar the surface density 
increases. Colour helps in distinguishing the main features 
supported by the orbits in our model. It is clear that, de- 
pending on the family which prevails, we could have the 
following cases: (a) If the ring consists of s, svl and sv3 or- 
bits at low energies it will be an oval with a more or less 
strong lemon shape, i.e. a morphology similar to the one 
of the NGC 6782 inner ring, (b) If the families of the f- 
tree are dominant, they will lead to the appearance of more 
hexagonal-like oval structures with corners close to the mi- 
nor axis (like in UGC 12646), (c) A combination of all oval- 
supporting families, leads to the lemon-shaped oval, together 
with segments parallel to the minor axis and corners close 
to the minor axis of the bar. This is the morphology of the 
inner ring of IC 4290 and is represented in most of its details 
in Fig. ITUt . 

For the rings it is crucial that their width (in the equa- 
torial plane) is not too large and this is indeed the case in 
Fig. 110b . In such a blurred image the filtering simulates the 
effect of trapping around the relevant periodic orbits. To 
show that the width of the gaussian that we have adopted 
for the filtering is reasonable, we examined the ring structure 
formed by single regular orbits, trapped around the periodic 
orbits discussed above. We found that it is well confined on 
the face-on view of the model, as Fig. Illl shows. This is due 
to the nature of the characteristic diagram in the region un- 
der consideration, since there a small deviation in the initial 
conditions is sufficient to stop the orbit from being trapped 
around the ring-producing periodic orbits. I.e. the trapping 
around the periodic orbits, whenever it happens, is of rather 
small extent. This means that the rings composed of such 
orbits will have small width, in good agreement with the 
observations. 

In Fig. 112b we give another face-on weighted profile con- 
sisting of stable orbits of the families f at energies —0.1990 < 
Ej < -0.1978 and fv5 at energies -0.1993 < Ej < -0.1969. 
In Fig. [T2"b we add also orbits from s for —0.1984 < Ej < 
-0.1972 and sv5 for -0.1984 < Ej < -0.1972. This combi- 
nation gives another type of rings, which can be described 
as polygonal. The orbits, as we have seen in the descrip- 
tion of the individual families, are of 6:1 and 8:1 types. The 
overall polygonal morphology, in Fig. 112b . however, can be 
described as hexagonal. The inner ring of NGC 3351 could 
be explained based on these families. 

The morphologies in Figs. [TJ;,d can also be explained 
with the stable orbits of the region, but only as specific 
cases, contrarily to the many possible combinations of or- 
bits we can use for building the inner ring morphologies in 
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Table 1. Energy intervals for the stable orbits contributing to the weighted profiles in Fig, 1101 and Fig. 1121 
The successive columns give the name of the family and the corresponding energy intervals in the figures: 



family 


10a 




10b 




10c 




12a 


12b 




s 


-0.1990 -» 


-0.1985 


-0.1990 -+ 


-0.1985 


-0.1990 -» 


-0.1985 




-0.1984 -4 


-0.1972 


svl 


-0.1989 -> 


-0.1968 


-0.1989 -» 


-0.1968 


-0.1989 -> 


-0.1968 








sv3 


-0.1987 -> 


-0.1982 


-0.1987 -> 


-0.1982 


-0.1987 -» 


-0.1982 










-0.1976 -» 


-0.1958 


-0.1976 -+ 


-0.1958 


-0.1976 -» 


-0.1958 








sv5 
















-0.1984 
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-0.2016 -+ 
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-0.2016 -> 


-0.1993 


-0.1990 -* -0.1978 


-0.1990 -4 


-0.1978 


fvl 






-0.2016 -» 


-0.1974 


-0.2016 -> 


-0.1974 








fv3 






-0.2002 -» 


-0.1967 


-0.2002 -> 


-0.1967 
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-0.1969 




I 11 I 1 11 





Figure 10. Weighted profiles obtained from stable orbits of families s, svl and sv3 (a) and s, svl, sv3, f, fvl and fv3 (b). In (c) we 
give a blurred representation of (b) and we see clearly a polygonal structure embedded in the ring. Colour in this last frame gives the 
projected surface density. The orbits are chosen at the appropriate energy intervals in order to support an oval shaped ring. 





Figure 12. Weighted profiles obtained from stable orbits of fam- 
ilies f and fv5 of the f-tree (a), and s, sv5, f, and fv5 (b). These 
are orbits at larger energies than the orbits of the s, f families 
used for the weighted profiles in Fig. 1101 
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Figure 11. A single non-periodic orbit trapped around the stable 
s periodic orbit at Ej = —0.0198176. It is integrated for more 
than 20 orbital periods. It forms an oval of small width, the inner 
boundary of which is lemon-shaped. 



Figs. ^,b. A 5:1 morphology (Figs. 0j) can be explained 
by combination of orbits belonging to frl (Fig. Qp), or to 
fvl.l (Fig. or to both. There are, however, two prob- 
lems. First, as we noted before, we can have a pentagonal 
symmetry when only one of the two symmetric branches of 
the family is populated, which presumably necessitates spe- 
cific initial conditions. The second is that these families have 
small stability regions. Both are bifurcated as stable and be- 
come unstable for an energy close to the bifurcating point. 



10 Patsis et al. 



For these two reasons pentagonal rings like in the case of 
NGC 3367 should be rare. 

Hexagonal rings with cusps on the major axis of the bar 
can be formed by material trapped around periodic orbits 
near the energy minimum of the tangent bifurcation which 
brings family s in the system. Nevertheless, stable s orbits 
with the morphology shown in Fig. [SJi exist only in a very 
narrow energy interval. The transition to the morphology of 
Fig- EH 3 happens very fast, while the orbit in Fig.||J is an ex- 
ample from the unstable branch of the tangent bifurcation, 
which remains unstable for a large energy interval. Indeed, 
these orbits become stable only when they have developed 
large loops. The ra rity of inner rings of the hexagonal type 
seen in NGC 7020 feutalll99oh reflects the way family s is 
introduced in the system through the tangent bifurcation. 
The stable branch of this bifurcation is associated with the 
formation of ovals, while the unstable one is associated with 
'NGC 7020-type' hexagons. The morphology of the stable 
and the unstable orbits is very similar only very close to the 
energy minimum of the s characteristic. Thus, the hexag- 
onal geometry we observe in NGC 7020 can be due only 
to non-periodic orbits trapped in this small fraction of the 
phase space. Note that the inner side of many 2D and 3D 
banana-like orbits trapped around the Lagrangian points 
L4 and L5, consider ed in pairs, can also form hexagons (see 
ISkokos et al"ll2002al - Figs. 18 and 19). They can not, how- 
ever, explain rings, because of their orientation. If we would 
consider at a given energy the two stable banana-like orbits 
trapped around L4 and L5, respectively, we build a ring elon- 
gated along th e minor axis of the bar rather than around 
the major axis JPatsis et alJI 20031 Fig.2), which is contrary 
to observations. 

Finally we note that NGC 6782 gives also an example 
of diamond morphology in the bar, due to orbits with that 
shape, which is not related to the ring. This can be realized 
by inspection of the innermost isophote drawn in Fig. ^1. 
These isophotes may correspond to xl diamond-shaped or- 
bits. 



6 CONCLUSIONS 

We have studied the orbital structure of inner rings in static 
3D Ferrers bars. The stable orbits we find belong to families 
which are typical for barred potentials with a type-2 gap 
in their characteristi c at the radial 4:1 resonance . In 2D 
jAthanassoulalll992aTl and 3D JSkokos et al.ll2002bl) Ferrers 
bars this is the usual type of gap we encounter. So we expect 
our conclusions to be valid for a large area of the parameter 
space. This work is the first step in a study of ring structures 
encountered in time-dependent, fully self-consistent iV-body, 
as well as in gaseous-response models. 

The main conclusions of the present study are the fol- 
lowing: 

(i) In 3D Ferrers bars, inner rings are due to orbits be- 
longing to families in the upper part of the type-2 gap at 
the inner radial 4:1 resonance. They are grouped in two or- 
bital trees, which have as mother-families the planar f and s 
orbits. For building the rings one cannot invoke orbits from 
the xl-tree or other families. The orbits that make the rings 
belong in their vast majority to three-dimensional families 
of periodic orbits. 



(ii) The 3D bifurcating families of the two groups (fvl, 
fv3, fv5; svl, sv3, sv5) play a crucial role in the morphology 
of the inner rings. They have large stable parts and thus 
they increase considerably the volume of the phase space 
occupied by ring-supporting orbits. The energy width over 
which we can find stable 3D orbits supporting the rings is 
larger than the corresponding interval of 2D stable families. 

(iii) The prevailing types of inner rings are variations of 
oval shapes and are determined by the way the f and s fam- 
ilies are introduced in the system, i.e. by the tangent bi- 
furcation mechanism. The orbits on the stable branch of 
their characteristic, together with their stable 3D bifurca- 
tions, support ovals with a more or less strong lemon shape, 
or oval-polygonal rings with 'corners' along the minor axis 
of the bar. These types of inner rings represent frequently 
observed morphologies. 

(iv) Pentagonal rings are rare because the families build- 
ing them (frl, fvl.l) have small stable parts and usually 
come in symmetric pairs. Thus, in order for these rings to 
appear, the symmetry must be broken and only one of the 
two branches be populated due to some particular forma- 
tion scenario. Furthermore, considerable material should be 
on regular non-periodic orbits trapped around stable peri- 
odic orbits existing only in narrow energy ranges. 

(v) If orbits are trapped around stable s periodic orbits 
at the energy minimum of the s characteristic, then an 
NGC 7020 morphology can be reproduced. Although such 
a morphology is in principle possible, it should be rare, be- 
cause it would necessitate that considerable amount of ma- 
terial be on regular orbits trapped around periodic orbits in 
a very narrow energy interval. Indeed the hexagonal orbits 
with cusps on the major axis are on the unstable branch of 
the tangent bifurcation. 
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